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Abstract 

It is known that a viscous fluid flow with curved streamlines can support 
both Tollmien-Schlichting and Taylor-Gortler instabilities. The question of 
which linear mode is dominant at finite values of the Reynolds numbers was 
discussed by Gibson and Cooke (1973). In a situation where both modes 
are possible on the basis of linear theory a nonlinear theory must be used to 
determine the effect of the interaction of the instabilities. The details of this 
interaction are of practical importance because of its possible catastrophic 
effects on mechanisms used for laminar flow control. Here this interaction is 
studied in the context of fully developed flows in curved channels. Apart from 
technical differences associated with boundary layer growth the structures 
of the instabilities in this flow are very similar to those in the practically 
more important external boundary layer situation. The interaction is shown 
to have two distinct phases depending on the size of the disturbances. At 
very low amplitudes two oblique Tollmien-Schlichting waves interact with a 
Gortler vortex in such a manner that the amplitudes become infinite at a finite 
time. This type of interaction is described by ordinary differential amplitude 
equations with quadratic nonlinearities. A stronger type of interaction occurs 
at larger disturbance amplitudes and leads to a much more complicated type 
of evolution equation. The solution of these equations now depends critically 
on the angle between the Tollmien-Schlichting wave and the Gortler vortex. 
Thus if this angle is greater than 41.6° this interaction again terminates in a 
singularity at a finite time; otherwise the breakdown is exponential taking an 
infinite time. Moreover the strong interaction can take place in the absence of 
curvature, in which case the Gortler vortex is entirely driven by the Tollmien- 
Schlichting waves. 

This work was supported by the National Aeronautics and Space Admin- 
istration under NASA contract No. NAS1-18107 while the authors were in 
residence at ICASE, NASA Langley Research Center, Hampton VA 23665. 
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1. Introduction 

In recent years there has been renewed interest in the development of 
laminar flow wings. In particular the wing developed at NASA Langley, see 
for example the paper by Harvey and Pride (1982), has been designed so 
as to achieve laminar flow by a combination of suction and passive mecha- 
nisms. This particular wing has significant regions of concave curvature on its 
lower side so the possibility of transition through the Gortler vortex instabil- 
ity mechanism exists. However, the complicated three-dimensional flow over 
the wing can also support Tollmien-Schlichting waves and crossflow vortices 
so the transition process for the flow is likely to be dominated by an inter- 
action involving all three instability mechanisms. Thus an understanding of 
the nonlinear interaction of these instabilities is necessary if they are to be 
successfully supressed. 

Here the nonlinear interaction of Tollmien-Schlichting waves and Taylor- 
Gortler vortices in fully developed flows in curved channels is discussed. Here- 
after we refer to these modes as “TS” waves and “TG” vortices respectively. In 
particular the interaction between lower branch TS waves and small wavenum- 
ber TG vortices is investigated. It is this type of interaction which will in 
general occu r in external flows, bu t with the added complication of boundary 
layer growth. It is known from the work of Hall (1983) that there are signifi- 
cant differences between TG vortices in internal and external flows; the major 
difference is of course that non-parallel effects cannot be ignored in the exter- 
nal case. Such a difference does not exist between TS waves in external and 
internal flows and the reader is referred to the papers by Smith (1979a, b), Hall 
and Smith (1982, 1984) for a discussion of the linear and weakly nonlinear 
stages of the growth of these instabilities. More recently Smith and Burggraf 
(1985) and Smith and Stewart (1987) have given a description of the nonlinear 
development of TS waves well beyond the weakly nonlinear stage. 

Though the primary motivation for the present calculation comes from 
the external boundary layer problem there are many engineering situations 
where fully developed flows in curved channels can support both types of insta- 
bility . In addition further motivation for a study of this interaction problem 
comes from the experimental results of Nerem, Seed and Wood (1972) who 
found that blood flow in the curved part of the aorta can support turbulent 
bursts for part of the cardiac cycle. It is known from the work of Sobev 
(1978) and Papageorgiou (1987) that TS and TG instabilities can occur in 
unsteady flows in curved pipe flows so it is likely that the transition process 
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there is again dominated by an interaction between these types of distur- 
bance. Experiments of some relevance here , in curved channels or boundary 
layers ,are presented by Peerhossaini , Clement and Wesfried (1985) ,Ko- 
hama(1987) ,and Winoto and Crane (1980) . 

The only previous investigation of the external interaction problem is 
that due to Nayfeh (1981). In that calculation Nayfeh arbitrarily assigned an 
amplitude to a Gortler vortex obtained by solving the parallel flow instability 
equations. The effect of this disturbance on a pair of linear oblique TS waves 
was then calculated. Nayfeh found that the vortex could have a large effect 
on the growth rate of the TS waves. However, more recently Malik (1986) 
repeated Nayfeh’s calculation and found that the effect of the TG mode on 
a TS wave was negligible. The difference was apparently due to an inconsis- 
tency between Nayfeh’s Orr-Sommerfeld boundary layer scaling and that for 
the vortex flow. Notwithstanding the doubts about the conclusions reached 
by Nayfeh on the basis of his numerical calculation there are other reasons 
why it is possible that his work is fatally flawed. The most obvious reason 
why Nayfeh’s approach is not valid is that it is known from the work of Hall 
( 1982a, b) and Hall (1983) that in the wavenumber regime appropriate to his 
calculation, a linear TG vortex is fully non-parallel and the parallel flow the- 
ory used by Nayfeh gives results having no connection with the non-parallel 
theory. Furthermore, the amplitude of the TG mode cannot be specified 
arbitrarily; it must be calculated by solving the Navier-Stokes equations nu- 
merically as in Hall (1987). The only situation where this is not necessary 
is when the amplitude of the TG mode is sufficiently small for linear theory 
to be applicable. However, in that situation the TS and TG modes become 
coupled so that again it is not possible to specify the amplitude of the Gortler 
vortex. 

In internal flows it is possible to study the two types of interaction prob- 
lem alluded to above. Firstly, it is reasonably straightforward to calculate 
numerically the fully nonlinear three-dimensional flow in a curved channel 
even at Gortler numbers much above the critical. The linear instability of 
this flow to linear TS waves can then be investigated in a self consistent man- 
ner. This type of calculation has already been carried out by Bennett and 
Hall (1987), who found that both two and three-dimensional TS waves can 
be greatly destabilized at relatively low Gortler numbers. Here the second 
type of interaction problem in which small TS and TG disturbances interact 
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with each other is studied. In neither of these approaches is it possible for 
the amplitude of the TG vortex to be specified arbitrarily. 

The interaction problem which occurs at very low levels of the distur- 
bance amplitude can in the first place be most conveniently studied through 
the nonlinear equations governing fully nonlinear unsteady TG vortices in 
fully developed channel flows. This is surprising since these equations do 
not contain the interactive pressure gradient in the flow direction, yet it is 
well-known that this gradient is crucial in determining the structure of two- 
dimensional TS waves. However the interactions which are considered in this 
paper involve oblique TS waves and there is sufficient structure in the non- 
linear TG vortex equations to enable unstable three-dimensional TS waves 
to develop also. The oblique modes initially considered have their wavefronts 
almost parallel to the flow direction. Later this restriction is relaxed and TS 
waves inclined at an 0(1) angle to the flow direction are investigated. 

The first weak interaction discussed above involves a TG vortex and two 
oblique TS waves inclined at an equal angle to the flow direction. The in- 
teraction is found to be characterised by three coupled quadratic amplitude 
equations typical of resonant triad interactions. In general it is found that 
this weak interaction leads to a finite time singularity in the solution of the 
amplitude equations. This type of interaction will be referred to as a type 
I interaction. The form of the solution when the singularity appears sug- 
gests a second type of stronger interaction, hereafter referred to as a type II 
interaction. In this situation the amplitudes of the TG and TS modes are 
much larger and the interaction is governed by a partial differential system 
for the TG mode coupled to an integro-differential equation for the TS wave. 
A major property of this type of interaction is that it can take place in a 
straight channel; in this case the longitudinal vortex is driven by the interact- 
ing oblique TS waves. In addition the longitudinal vortex is sufficiently strong 
so as to react back on the TS waves. Here again the ultimate effect of the 
interaction is to cause the growth of both disturbances to a more nonlinear 
state. However the nature of this breakdown now depends crucially on the 
orientation of the oblique TS waves. In particular it is shown that TS waves 
propagating at an angle of more than 41.6° to the flow direction produce a 
breakdown at a finite time. TS waves inclined at an angle of less than 41.6 
degrees also lead to a breakdown but this takes place over an infinite time 
interval and is therefore less catastrophic. Each type of breakdown for these 
type II interactions can occur with or without curvature of the channel; in 
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the presence of curvature the route to breakdown can be more complicated 
but the final state has little dependence on the curvature. 

The procedure adopted in the rest of this paper is as follows . In §2 the 
nonlinear equations governing the growth of nonlinear TG vortices in fully 
developed channel flows are formulated. In §3 these equations are solved for 
weak type I interactions, whilst in §4 the corresponding equations governing 
much stronger type II interactions are derived. In §5 the solution of the in- 
teraction equations found in §4 is discussed. Finally in §6 the results obtained 
in the previous sections are discussed and some conclusions drawn. 

2. Formulation of the nonlinear disturbance equations 

Consider the steady flow of a viscous incompressible fluid of density p*, 
kinematic viscosity v* in a curved channel driven by an azimuthal pressure 
gradient. The channel walls are defined by r' — R\, r' — R2 with respect to 
cylindrical polar co-ordinates ( r',0',z ') and it is assumed that the curvature 
parameter 8 defined by 

5 = ( 2 . 1 ) 

is small. The co-ordinates [x, y, z) are defined by 

RiO' _ r'-R 1 _ z' 

x ~ R2-R1’ v ~ R 2 - Ri ’ z ~ R 2 -Ri 

and the corresponding velocity field ( u,v,w ) in the small gap limit 8 « 1 
can be written 


(u,v,w) = u B = t/o(«(y),0,0) +O(8U 0 ) (2.2) 


where 

u = y{ l-y), 


whilst is the maximum value of the zeroth order flow velocity. The insta- 


bility of (2.2) to three- dimensional time-dependent Taylor-Gortler vortices 
and Tollmien- Schlichting waves will now be considered. 


The appropriate non-dimensional time scale t for the disturbances is slow 
scaled on (R 2 — Ri)/Uo, whilst the perturbation velocities have the usual TG 
scales. The total velocity field is now written as 
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U = 11 B + UoiU^E^R^W), 
where Re is a Reynolds number defined by 


Re = 


Uo(Ri - Ri) 


(2.3) 


(2.4) 


and U,V,W and the pressure perturbation P, scaled on p*u* 2 (R 2 — Ri) 2 , 
are all functions of x, y, z and t. The Gortler number G is defined by 


G = 2R 2 e 6 


(2.5) 


and is taken to be 0(1) in the limit <5 — >• 0. The disturbance equations correct 
to 0(6°) can then be written in the form 


£(u, v, w) = (0, P y , P z ) + (u.V)u + (vu\ Guu + G^-, 0) 


Ux “ 1 “ Vy ©2 0 ) 


(2.6a, 6, c,d) 


where 


r = — — - 

~ dy 2 dz 2 dt 


It should be noted that the pressure does not enter the x-component of the 
momentum equations; this is the most significant difference between the equa- 
tions for TG modes and TS modes inclined at an 0(1) angle to the basic flow 
direction. However, the equations (2.6) have sufficient structure to support 
TS modes which have wavefronts inclined at a small angle to the flow direc- 
tion (see Hall and Smith (1984)); it is the interaction of these modes and the 
TG mode which will be first investigated. 

The equations (2.6) reduce in the small amplitude case with = 0 to 
Dean’s (1928) equations governing the linear instability of Poiseiulle flow in 
a curved channel. It is known that instability occurs for sufficiently large G 
and that (3 , the z - wavenumber of the vortices, is proportional to either G~ * 
or G* where G » 1. The manner in which nonlinear effects operate in these 
limits is quite different, thus for (3 « 1 the interactions are governed by a 
Stuart- Watson type of interaction whilst for (3 » 1 a mean field interaction 
of the type discussed for external flows by Hall (1982b) applies. Here the most 
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relevant case is when (3 « 1 so that there is a strong interaction between the 
TG and TS modes. The nature of the interaction depends on the size of the 
disturbance velocities; in the next section the simplest nonlinear interaction 
of the resonant triad type is considered. Later more complicated higher- 
amplitude interactions which lead to partial differential evolution equations 
are considered. 

3. Type I interactions between Tollmien-Schlichting waves and 
Taylor- Gortler vortices. 

Suppose that a TG mode of instability with cross stream wavenumber 
/ 3 is superimposed on the basic state (2.2). If (3 is small the disturbance 
equations (2.6) can support TS waves inclined at a real angle 0(/? 7 ) to the 
flow direction. Thus the pair of skewed TS waved proportional to 


Ef- 1 = exp 


E^ 1 = exp ±z 


ax 


(3 6 


i/3z 

~Y 


ax i(3z 


will interact nonlinearly with a TG vortex proportional to 


Ef 1 = exp ( ±i(3z ) 


(3.1a, 6) 


(3.2) 


at second order, to reinforce themselves and the TG mode. This resonant 
triad type of interaction has been investigated elsewhere (see, for example, 
Craik (1971) or Smith (1987)) in the context of purely TS instabilities of 
boundary layers. 

The crucial size of the disturbance associated with (3.1) and (3.2) is 
fixed from the following considerations. Firstly suppose that the TS and TG 
modes have amplitude A((3 ) and C{(3) for the streamwise velocity components 
in the core. The relative size of A and C is chosen such that the amplitude 
equation for A driven by terms proportional to AC arises at the same order 
as the equation for C driven by terms proportional to \A \ 2 so that nonlinear 
interactions occur over comparable time scales. Further limits away from this 
regime can then be taken to examine the linearized problems associated with 
A/C — > 0 or C/A — * 0. Before writing down the scalings deduced from the 
above considerations it is convenient to define a slow time scale r by 


r = (3 k ~ A t 
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and expand the Gortler number G in the form 

G = -j^{go + (3 k ~ 4 gi H — } 


(3.3) 


where k > 4 is a constant associated with the elevation of the Gortler number 
G above the neutral value go(3~ 4 . The expansion procedure given in this 
Section fails when k = 4, in which case the appropriate evolution equations are 
partial differential equations but can be derived by some minor modifications 
to the present procedure. It is now necessary to be more precise about the 
downstream wavenumber of the TS waves; thus E\ and E 2 are now taken to 
be 

„ . f ax idz i~it 

«i=«P*|^ + T 

i/3z i~ft 1 

) 

where 


En — 


= exp i | 


ax 

w 


I 


(3.4a, 6) 


= a 0 + + Qi 2P z h + 


^3.5a,6) 


7 = 7o +7i/? + 


The ordering of the term (3 k in (3.5a) depends on the size of k; we assume 
that all the quantities in (3.5a,b) are real so that 


jn ip __ ii 1 77* fj 1 §j 1 jp jp _ , _ TP 

Jz/lH/2 — il/ 3 , Hj 2^3 — ^1) JZ/lXL/3 — XI/ 2 * 


In the absence of a Gortler vortex (3.5a,b) would of course reduce to the lower 
branch asymptotic expansions for the TS neutral wavenumber and frequency 
at a given value of (3. 


The coreflow solution 

The above considerations about disturbance size lead us to write down 
the following expansions for the disturbance velocity field in the main part of 
the channel: 

U = fd k u 3 E 3 + + /3 k+3 u 2 E 2 + • • • + C.C., 

V = p k v 3 E 3 + /3 fc " 3 t;i£i + p k ~ 3 v 2 E 2 + • • • + C.C., 

W = pb-'wsEs + /3 fc “V-Ei + P k ~ 2 w 2 E 2 + • • • + C.C., 

P = fi k ~ 2 p 3 E 3 + /? fc - 9 pi£i + f3 k ~ 9 p 2 E 2 + • • • + C.C. 


(3.6a, 6, c, d ) 
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where 

u i = u io + **• + /? uu H » 

U 2 = U 20 + • • • + 0 k U2i + • • • , (3.7a, 6, c) 

= tig 0 -|- (3 k usi + • • • . 

together with similar expansions for Vi , V 2 , V3 , w 1 , W 2 , «>3 , pi , P 2 , P3 . The terms 
in (3.7a,b) represented by • • • depend on the particular size chosen for k but 
do not enter the calculation given in this Section. The TS modes in (3.6) are 
valid only away from viscous layers of thickness (3 2 at y = 0, 1 where the no 
slip condition must be satisfied. In contrast the Gortler expansions are such 
that viscosity enters the coreflow problem at zeroth order and hence the no- 
slip condition can be satisfied directly. The disturbance expansions (3.6), (3. 7) 
are then substituted into (2.6) and like powers of (3 are equated. The zeroth 
order problems for ( Uio,vio,wio,pio ) and (V30, «>30>P3o) are found to be 

ia 0 uu 10 + v 10 u' = 0, 

UkoUVio + PlOy = 0, 

. _ i 

lOtOUWio = — -P10, 
icxQU 10 + Vioy = 0. 

for the TS part and 

u 30yy — V30U y , 

P 30 y = 90UU30, 

^30yy — ®P30> 

v 30y + iwzo = 0 . 

for the Gortler part. 

As expected the system (3.8) for the TS part has no solution with the 
no-slip condition at y = 0, 1 is satisfied. Thus if the boundary conditions on 
uio , w 10 at y = 0, 1 are dropped the appropriate solution is 


(3.8) 


(3.9) 


P10 = 


«io = ao(r)u r , 



v\o = —iotoao(r)u, 


W 10 


~P 10 
2 a 0 u 


(3.10) 
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where a(r) is an as yet unknown amplitude function. A similar analysis shows 
that 

u 2 q = & 0 (r)u', v 2 o = ~ia 0 b 0 (T)u, 



(3.11) 


where bo is the second TS ‘amplitude’, or ‘displacement’ function to be found 
at higher order. 


The zeroth order Gortler equations (3.9) retain the viscous derivatives 
in the y direction and therefore can be solved subject to u 3 q = v 30 = W 30 = 
0, y — 0, 1. However this system will have a solution for only certain values 
of go] the lowest one of which corresponds to the left hand branch of the 
lowest neutral curve for the Dean problem. The solution can be written 


^30 = co (r) C/30, v 3 o = co(r)T 3 o, W30 = co(t)Wso, P30 = co(r)P 3 o 

(3.12a, 6) 

where c 0 (r) is another amplitude function and U 3 0,^30 are determined by the 
eigenvalue problem 

U30 = V30U 1 , 

V 30 ~ 90 UU 3 Q — 0, (3.13) 

C/30 = V ' 0 = V 30 = 0, y = 0,1. 

The functions un, u 2 i, U31 etc. are found by proceeding to higher order in the 
coreflow expansions. In particular it is found that (tfcn, vu,wupu) satisfies 
an inhomogeneous form of (3.8) forced by the interaction of the TS waves 
proportional to E 2 and the Gortler vortex. The appropriate solution of the 
system is 

uu = a!(r)u' + 6 0 c 0 CC'o 

Vn = —iotoa\{r)u — ic*o&oCoC/3o 
ry rv 

P 11 = -aoai(r) / u 2 dy-2alb 0 c 0 / U 3 oudy, 

Jo Jo 

where ai(r) is another amplitude function to be determined at higher order. 
The solution for ( u 22 ,v 22 ,p 22 ) is given by (3.1) with bo replaced by c 0 and a± 
by the function 61 (r). Finally to complete the coreflow solution to this order 
the functions u 3 i,v 31 ,u; 3 i and p 3i satisfy 


u 31yy — v 31 u ’ = «30 t ? 
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0 = ~P3iy + 9 ouu 3 i + guu 30 - ^a 0 b 0 otQUu\ 

w 31 yy — *P 31 + ^ 30 t + ia^QO^U 2 , 

^ 31 y + iw 31 = 0. 

Here the forcing terms arise from the interaction of the two skewed TS waves 
in the core. An examination of the TS waves in the wall layers shows that the 
interactions in the wall layers have a negligible effect on the Gortler vortex in 
these regions so that the above sysytem is to be solved subject to the no-slip 
condition at y = 0, 1. In this case a suitable solution exists if and only if an 
appropriate orthogonality condition is satisfied; this condition can be written 

^ = g$c o + ~ia 0 b 0 (314) 


where 

g _ ~ Jo 1 V3Q*P(y) dy 

fo(U so Q(y)+V&P(y))dy’ 

and 

_ 2q:q Jq 1 uu'P(y) dy 

1 ~ ti(V 3 oQ(y)+V;' 0 P(y))dy' 

Here P,Q are solutions of the adjoint problem 

Q” ~ 9ouP = 0 

P"" = Qv! 

with P = Q = P' = 0, y = 0, 1. 

At this stage in the coreflow problem the TS amplitude functions are not 
determined and in fact are found by solving the equations in the wall layers. 

The solution in the wall layers. 

At the upper and lower boundaries there exist layers of thickness f3 2 where 
the TS waves adjust their structure so as to satisfy the no-slip condition. Thus 
near y = 0 it is necessary to define a stretched variable Y by 


Y = 0~ 2 y 
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and (3.6), (3.7) are modified by writing 


u = /3 2+ *y E 3 c o U^(0) + + U 2 E 2 } + • • • + c.c. 

V = 0*+ k Y 2 E 3 c, o V 3 "(0) + + V 2 E 2 } + ■ ■ ■ + C.C. 

w = 0 k+1 YE 3 coW^O) + ^-'{WiE! + W 2 E 2 } + • • • + C.C. 
p = 0 k E 3 c o P 3 (O) + 0 k ~ 9 {P 1 E 1 + P 2 E 2 ] + . . . + C.C. 


where 

Ui = U 10 H 1- /3 k Un + • • • , 

V 2 = U 20 + • • • + (3 k U2i + • • • . 

Here the first terms in the expansions (3.15a,b,c,d) are just the first terms 
in the small y expansions of the zeroth order Gortler vortex in the core. 
Again the ordering of the correction terms of order (3 k in (3.16a,b) depends 
on the precise value of k. For brevity the wall layer solution for the TS wave 
proportional to E\ will be given, the details for E 2 being identical apart from 
s o me si gn and n otational chang es. The expaii3ions_(3J.5),_(3Ji)) are th en 
substituted into (2.6) and like powers of (3 are equated. At zeroth order it is 
found that (Uio,Viq,Wio,Piq) satisfies 


(3.16o,6) 


{-^-i<*oY + i-io}V 10 = V lo 

dP 10 


dY 


= 0 


d 2 


(3.17a, 6, c, d) 


- ia 0 Y + i l0 }W 10 = 2^10 

iotoUio + Vio + — Wio — 0 

and these must be solved subject to the conditions U 10 = Vio = Wio = 
0, Y = 0. 

The appropriate solution is obtained in terms of the variable £ defined 

£ = A {Y-H}, A = (ia 0 )* 

OiO 


by 
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and is given by 

a 0 U 10 + i/2W 10 = A 0 I Ai(rf) dr), 

* / £° (3.18a, 6) 

W l0 = lp 10 A~ 2 £(£). 

Here A, is the Airy function whilst £(f) satisfies an inhomogeneous Airy 
equation with unity on the right hand side and zero boundary conditions at 
0, oo. 

The amplitude function Aq can be related to the pressure Pio using the 
x and z momentum equations evaluated at Y = 0 to give 

= (ia a )i A[(Zo)Ao- (3.19) 

The zeroth order eigenrelation for the Tollmien-Schlichting wave follows by 
matching (3.18) and the corresponding upper layer solution with the coreflow. 
If we denote the functions in the upper layer corresponding to Uio,Aq etc. 
by U\q,Aq etc. it follows from (3.11) that 


Pio — Pio 

a o r A , 

a 0 = — / Ai dy 

a ° Jz 0 

and using (3.19) and the corresponding upper layer relationship it is found 
that 

A 2 4(e°) = ^ r A ' Mdri (321) 

J £o 

and o:o,7o must be chosen to satisfy the above eigenrelationship. 

A similar analysis applies to the second TS mode which has amplitude 
Bq and clearly this mode leads to the same eigenrelation (3.21). The ampli- 
tude equations for ao(r) and &o(r ) will now be obtained by considering the 
0(/3 k ) disturbance equations in the wall layers. The equations to determine 
{UiuVn,Wu, P n ) in the layer at y = 0 are: 

d 2 

{-jyt-icxoY+^Uu-Vu = U^O)Yco{U 2 oiaoAiW 2 o}+V 2 oCoU^O)+iaYU 1 o+U 1 o 


<*6 „ (+\ 

~ “57r°o(*)» 


30 

—Aq ^ 


ao 


j 

Ji 0 


(3.20a, 6) 


Ai dr) 
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*2 

{jjTj - ia 0 Y + i l0 }W n - lp„ = C/'(0)yc o lV2o + iaYW w + W i0 r 
iaoVu + V 11Y + Wu = —iocUio 

Uu = Vn = W n = 0 ,y = 0 . (3.22) 

The appropriate solution of (3.22) is given by 

«*o{aptfu + \Wu} = »aoylo{i^.-(0 - fo^.({)} 

dAn ft 

+ A-^{A,(£) — A,(£o)} + »aoAi / Ai d£ 

J £o 

- coBb^(0)tooil5(fo){£(0 + ({- 2f 0 )£'(fl + &£'(&)} 
+ (/UO)m 0 coB 0 {^A,U) - f 0 M() + | (oA,((o)} (3.23) 


iP 
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A _ l{ 2i|Ao[ A ,( & ) _ + A^&Ai(f 0 ) 

+ iaoA l A[(£o) H — c o Po^3(0)[— lo^i(^o) + -4((£o)] 


- A'^coiaoBoU^O) 


3 + 


egx 

MZo)\ 


} 


(3.24) 


rOO 

Here A\ is another amplitude function to be found and X = Ai(rj)drj. 

Jto 

A similar analysis applies in the upper wall layer and the matching between 
the velocity and pressure fields where the wall layers merge with the coreflow 
produces an equation to determine ao(r) the amplitude of the first TS mode 
in the core. The equation, after some manipulation, can be written 


where 


^ - = aua 0 + nb 0 c 0 , 
dr 


V = r t 


2A ^(&) [^fey + i 


(3.25) 
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and 


H ={Aia 0 Jx + iaoK(O) - tf£(l)]{I [|{ 0 ili(&) 


+ &^X]_ 1 _|( 1 _ € g^|W) 

Ai(fo) -I 3\ AJ(f, o)/ 

+ ( 3+ ^))} 2AAi(?o) [lRlo) + ^] } 


with 


J = [ Uzo{y)u{y)dy . (3.26a, 6, c,) 

Jo 

A similar analysis for the second TS mode produces the amplitude equation 


d&o 

dr 


= ai/6 0 + Ma 0 c 0 . 


(3.27) 


Thus ao,6o> and Co satisfy the resonant triad equations (3.14), (3.25), (3.27). 
These equations can be put in a more convenient form by writing 


cio 



t 


= ah, 

•O 

II 

o 

pO 

Co — cl 

/ry | e'*, 

k=l 

\h | 2 
\«9 \ ’ 

II 

T 

1^1 


v=\v\ e tw , (3.28) 


to give the equations governing the nonlinear interaction of TS and TG in- 
stability modes: 


da 

dt 

db 

dt 

dc 

dt 


ae iU a + e'^bc 
ae iU b + e'^ac 
±c + ab. 


(3.29a, 6, c) 


Here the ± signs in (3.29c) are taken depending on whether g > 0 or g < 0. 
The constants u and <f> can be found by evaluating the coefficients in (3.14), 
(3.25), (3.27) numerically and using (3.20). It was found that 


e i<f> = (0.998,0.049) 


(3.30a) 
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corresponds to go = 4.27 x 10 5 , the smallest neutral value of the scaled Gortler 
number. The second mode for the Gortler problem has g 0 = 4.94 X 10 7 and 
in that situation 

e*^ = (0.79,0.61). (3.306) 

The equations (3.29) permit the particular solution a — b with c real and the 
following discussion is restricted to this situation. In this case it is convenient 
to write p =| a | 2 so that (3.29) can be replaced by 


~ {“ cos w + c cos 
dc 

Tt= ±c + p - 


(3.31) 


Apart from the trivial solution p — c — 0 an equilibrium solution of (3.31) 
corresponding to a stationary Gortler vortex and a time-periodic TS wave is 


p = rp c = ± 


a cos u 
c os 


(3.32) 


It will be seen shortly that (3.31) apply to another interaction problem in- 
volving TS waves and TG vortices so it is convenient to discuss (3.31) in 
general and later discuss the special case corresponding to (3.30). However in 
all the situations which arise in this paper cos u > 0 so that the properties of 
the limit cycle solution depend only on the sign of cos <ft. A linear instability 
analysis of the limit cycle solution leads to the following conclusions: 

Case a cos <ft > 0 

Here the limit cycle exists if the TS and TG modes are both linearly stable 
or both linearly unstable. In either situation the limit cycle is unstable. 

Case b cos <ft < 0 

Here the limit cycle exists if the TS mode is unstable and the TG mode 
is stable or vice versa. The former situation is stable and the latter situation 
is unstable. 

Another property of (3.31) which has some relevance is the question of 
whether the system has a solution which breaks down at a finite value of the 
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time t = t. In order to investigate this possibility p and c are expanded in the 
form 

p ~ it-t) 2 + ’ 

1 co ( 3 - 33a ’ 6 ) 

c - (* - 1) + ’ 

and substitution of these expansions into (3.31) yields 


1 = —Co cos </> 

1 

Po — —Co — T 

COS (f) 

so that a finite time breakdown occurs if cos <f> > 0. 

The discussion leading to the amplitude equations (3.31) is easily modi- 
fied to account for a related resonant triad interaction problem. It should be 
recalled that lower branch Tollmien-Schlichting waves inclined at an 0(1) an- 
gle to the flow direction have x and z wavenumbers scaled on R^ 1 ^ 7 . Since the 
z variation in Section 2 was taken to be on a length scale 0(1) it follows that 
the above situation will occur when = 0{R E 1 ^ 7 ) so that a = 0{R E 6 ^ 7 ). 
The expansion procedure leading to (3.31) is essentially unchanged in this 
region, the only significant alteration being that the pressure gradient in the 
x-momentum equation is no longer negligible in the description of the two 
skewed TS waves. The amplitude equation (3.14) for the TG mode remains 
unchanged whilst the nonlinear term in (3.25) has coefficient p given by 


2A(f 0 +^M) J 4,( & )^ = 


2c*0 J . -^t(£o) /rW /p>\ TTt ( l\\ 

- ia o -^-(U 3 {0) - U 3 (l)) 


A 2 


+ 


+ 


mo) ~ m]) - SHy + w 4i(&) ) ( “° ■ + t> 


2(8 2A'((o) , 2(o 


34(6>) . f 2X4(6>). 


' A(£o 


+ « 


+ 


f 4(eo)\$i 

o))t\ 


i(&) ' AUto) 
where (3q = flR E 1 ^ 7 and the eigenrelation corresponding to (3.21) is 


(3.34) 


A 2 4(&) = ^ + f)x 


(3.35) 


(NO 
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The amplitude equations can again be rescaled into the form(3.31) and 
the coefficients cos oj and cos <f> are now functions of (3 0 with c*o ? 7o determined 
by the eigenrelation (3.35). In Figure 1 the dependence of ao on fio is shown. 
For large /? 0 , ao asymptotes to zero like /3q 6 which is consistent with the 
earlier analysis of §3. The constant cos u is positive in all cases whilst e*^ is 
shown in Figures 2,3 for the first and second Gortler modes. The lowest mode 
always has cos <f) > 0 so that the finite time breakdown discussed above will 
always occur. The second mode has cos co < 0 and cos<f> < 0 for fio/ao < 1.29 
so that a stable finite amplitude equilibrium state is possible outside this 
region. However the first Gortler mode is unstable in this situation and might 
possibly destabilize this state. 

4. Larger amplitude Type II interactions governed by an integro- 
differential system. 

It is now shown that by increasing the disturbance size the coupled triad 
interaction equations derived in the previous section alter their form. The 
choice of k in (3.3) effectively fixes the size of the TS and TG modes in the 
core^ If fc > 4 t he zeroth order TG mode has as its zeroth order pr oblem ( 3.9). 
At this stage the interaction of the two TS modes in the core does not force 
the TG mode. However, if k is now taken to be 4 then this is no longer the 
case and the forcing terms which appear in (3.13) now appear on the right 
hand sides of the zeroth order problem corresponding to (3.9). In contrast 
the forcing of the TS modes by their interaction with the TG mode remains a 
higher order effect and does not enter the zeroth order eigenrelation problem 
involving 70 and ao- Thus if the analysis of the previous section is repeated 
with k — 4 throughout then, after eliminating w 30 and P30, it is found that 
U30 and V30 satisfy 


( d 2 d. d 2 

*<v <V V3 ° 

dvzo 

U30 = V30 = 


V30U', 


gouu 30 - 2 puu', 

0, y = 0,1. 


(4.1) 


Here attention has been restricted to the situation where the skewed TS modes 
have identical amplitudes and again p =| ao | 2 . Thus (4.1) now replaces (3.9) 
and (3.13) and the outcome is that the zeroth order Gortler problem is now 
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forced by the TS modes so that g 0 is no longer required to be an eigenvalue 
of the Dean problem. Indeed the structure that is now developing occurs in 
the absence of curvature in which case the streamwise vortices present in the 
flow are a direct result of the interaction of the TS modes. 

The derivation of the amplitude equation corresponding to (3.25) is es- 
sentially identical to that given in Section 3. However the appearance of U 3 
in (3.26) means that the corresponding amplitude equation for p =J ao | 2 now 
depends on w 3 o. In fact the appropriate equation is 


1 d /*1 

2 ^P = p{« l/ r + 7 / uu 3O dy + d{u 3O y(0,t)-u 3 oy{l,t)}} 

V 0 

where u is defined by (3.26a) whilst 

7 = Real (i fltoX /(2A.(fo)[-^ + i])}, 


(4.2) 


(4.3a, 6 ) 


9 = Real 


' *»[£(§&*(fe)) + - 1 - 1(1 - gjfjgg) + (3 + MjI 

2AA t -(f 0 )[£oMJ(fo) + xl 


If (3 is reduced to be 0(J? -1 / 7 ) then the three dimensional effects again 
modify the co-efficients in the amplitude equation (4.2) for the TS mode. For 
the purposes of this paper it is sufficient to note that when a = aoR ~ 6 / 7 , 
/ 3 = A )# - 1/7 the eigenrelation (3.35) applies but with 7 , 6 now defined by 


7 = Real 


| -agx(o§ + jf 

1 A3(f 0 + ^^)^(e 0 ) 


(4.4a, 6 ) 


6 = Real 


[{ 


-iotpA'^ &) ia 0 Ai ( £ 0 ) 


f (a 2 + 32 14] (i, 2 _ 

X (ag + /? 0 2 /4) X V ( ° + ^ o/ 3A t (6) 

1 ( 3 f\ ( H I 


2A;.(e 0 ) + 2e 0 A;(6) 


v^o)'^Hfor + w)4)}/ 2A( ^ + 


3 3 X 


) 


X 


) 
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For convenience (4.2) is rescaled by writing a = 2av r , 0 = 2 6, 7 = 27 
so that p now satisfies 

dp f 1 

— = p{a + ~ij uu 3 ody + 9[u3Oy{0,t)-U3Oy(l,t)]} (4.5) 

which must be solved together with (4.1). It should be stressed at this stage 
that setting go = 0 in (4.1) does not force U30 = V30 = 0 so that in the absence 
of curvature the longitudinal vortex is driven by the interaction of the skewed 
TS waves. 

The solution of the coupled system (4.1), (4.5) is discussed in the follow- 
ing section. It will be seen that the nature of the solution depends crucially 
on the ratio 6 / 7 which must be calculated numerically from (4.3) or (4.4b) 
depending on the wavenumber under investigation. In Figures 4a, b the values 
of 7 and 0/7 calculated from (4.4) are shown as functions of Po/oto- The ratio 
0/7 is seen to increase monotonically from —.0138 at Pol oto — 0 to .028 when 
Po/oto —*■ 00. This limiting value corresponds to the value of 0/7 which is 
obtained from (4.3). Th e dis cussion of the fo llowing section shows t hat th ere 
is a significant difference in the solution of (4.1), (4.5) from 0/ 7 < —1/140 and 
0/7 > —1/140. For future reference it should be noted that 0/ 7 = —1/140 
when Po/cto = .8887, so that there is a change in the type of nonlinear inter- 
action when the angle between the TG vortices and the TS waves falls below 
41.6°. 

5. Solution of the evolution equations for Type II waves. 

Initially it is assumed in the following discussion that 0 and 7 in (4.5) 
are arbitrary real constants; later solutions of (4.1), (4.5) for values of 0 and 
7 appropriate to TS waves in channels will be considered. Firstly it is noted 
that an equilibrium solution of (4.1), (4.5) is 

p = U 3 = v 3 = 0, (5.1) 

but an instability analysis of this equation shows that this solution is unstable 
if go or a is positive. In addition a limit cycle solution exists with ao periodic 
in time and p , U 3 and V 3 stationary. This solution takes the form 


u 3 = U(y), v 3 = V(y), p = p c 


(5.2) 
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where 


U" = Vu', V"" = —2p e uu' + U = V = V' = 0, y = 0, 1, (5.3) 


and 

a + 7 f uUdy + 0[U'(O)-U'(l)] = O (5.4) 

Jo 

For a given value of <70 it follows from (5.4) that, since p e is necessarily 
positive, the limit cycle solution exists for either a < 0 or a > 0. In fact 
the choice of sign required for the existence of the limit cycle solution is 
determined by 6/ 7. A numerical solution of (5.3) for different values of <70 
shows that ^ 

7 j uUdy + 0[U'(O)-U'(l)] = O (5.5) 

on the discontinuous curve C shown in Figure 5. The integral J = uU dy 
and the shear difference fi = U (0) — U* (1) change sign along C ; the different 
combinations of signs possible for 0 < <70 < 10 9 are indicated in that Figure. 
The Gortler numbers <701 and <702 correspond to the first two neutral curves 
at low wavenumbers in the absence of TS waves. At these positions /z and 
J necessarily change signs but these quantities are found to also change sign 
independently at different values of <70- These sign changes mean that the limit 
cycle exists for positive values of a when (0/7, go) lies below the continuous 
parts and above the dotted parts of C. Elsewhere the limit cycle solution 
exists for negative values of a. The stability properties of the limit cycle 
solution can be found from a small perturbation 


«3 = U(y)+ u(y)e 3t , v 3 = V(y) + v(y)e 3t , p— p e + pe 3t 
in which case the eigenvalue s is determined by 

( d* \ 


-/ 


dy 2 


— S ) V = vu , 


(ft \ (ft I) 

-s\^=g 0 uu-2puu\ 


\dy 2 J dy 2 

P s= J uudy + 0{u'(O) - «'(!)}) Pe 


(5.6) 
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It can be shown from (5.6) by expanding for | a \« 1 that the limit 
cycle solution is unstable whenever it exists for a < 0, | a |<< 1. Numerical 
investigations of (5.6) confirmed this result for a in this region and showed 
that the limit cycle remains stable for 0(1) negative values of a. These 
investigations also showed that the limit cycle solution is unstable whenever 
go > 0oi This is entirely consistent with the results obtained for the ordinary 
differential equations discussed in §3. This is to be expected since these 
equations can be retrieved from (4.1), (4.5) by taking the limit go — *■ goi . 

Thus the limit cycle solution is only stable for go < goi and values 9/ 7 
below C in Figure 5. For ordinary differential equations it is known that 
unstable limit cycles often lead to a finite time singularity beyond which the 
solution cannot be continued. Though the evolution equations (4.1), (4.5) are 
not ordinary differential equations the possibility of such a ‘blow-up’ is now 
investigated. 

In fact the nonlinear partial differential system (4.1)-(4.5) leads to a 
str onger singulari ty than t hat found in th e ordinary differ e ntial lo wer ampli- 
tude system of §3. It follows from (4.1) and (4.5) that if a singularity develops 
then the dominant balances in these equations have 


du? 


— v 3 , 


dv? 


p 


dp 

dt ~ “ 3P 


dt “ r - dt ’ 

which suggests that in the neighbourhood of the singularity at t = to,u 3 ,v 3 
and p expand as 


__ w 30 

u 3 — — re + 


Substitution of the above expansions into (4.1), (4.5) shows that 


V30 


P = 


P 0 


+ 


(5.7a, 6, c) 


U30 = -V 3 0 u', V30 = pouu', 

and the solution of the above system which satisfies u 3 = v 30 = 0 at y = 0, 1 
is 


P 0 

*“=T 


[/; 


W30 — 


‘dy-y J u 2 dy,j 

= ~irPo f u 2 dy-yfu 2 dy. (5.8a, 6) 

^ [Jo Jo 
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where u = y( 1 — j/). 

Thus it is not possible to satisfy v' 30 = 0 at the boundaries at this stage. 
This is of course a consequence of the fact that at 0(1) values of y the second 
derivative in (4.1) is always negligible compared to The remedy is to 
allow for extra sub-layers at y = 0, 1 and, as is usually the case for the heat 
equation operator Jp- — ^ , the layers are of thickness (to — t)%. Thus near 
y = 0 the variable £ is defined by 



y 

t to-t ]* 


(5.9) 


and (5.7) is then modified by writing 


u 3 - 


MO 

to — t I 1 / 2 


Vs = 


t 0 - t | S /2 


The function V 30 now satisfies 


P ~ 


P 0 


+ 


t 0 -t 

(5.10a, 6 , c,) 



td_ 

2 d£ 



= 0, 


and the appropriate boundary conditions are 


(5.11) 


V 30 = VL> = 0 , £ = 0 . 


The condition that (5.10b) matches with (5.8a) when £ — > 00 requires that 

^' 0 — fij. 

so that the appropriate solution of (5.11) is 


60^30 
P 0 


= f-vw(§,^) 


e«’/» 

U'( |, 0 ) + 


C/(|,0)s/2 

U'( |, 0 ) ■ 


Here U (a, x) denotes the parabolic cylinder function which decays when x — ► 
00 . The function £/ 3 o(£) then satisfies 


60 / rf 2 £ d 11 _ f _ U{ |,0)\/2 

Po l<2£ 2 2 d£ 2P 30- £/'(§, 0) U'( |,0) 
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Uao = 0, e = 0, 

whilst matching with (5.8b) requires that 


C7J 0 -► 

30 60’ 


oo. 


The appropriate solution for C/30 is 

60 e« s /8 2v ^ 1 r’ 5 _, rr ,l .. 


+ e+ 


2 ^(|, 0 ) 
U'(l, 0 ) 






(5.12) 


After some manipulation (5.12) yields 


^o(0) = _ 

so that the gradient in tt 30 at the edge of the boundary layer is exactly twice 
the value at the wall. An identical analysis to the above applies at the upper 
wall y = 1 and then (4.5) yields at zeroth order 


3 

2 



uu 30 dy + 20U^ O 


so that po is determined by 


3 _ Po [ 7 , .1 

2 60 L140 ' J 


(5.13) 


Therefore this type of finite time singularity can only occur whenever 


9 1 

7 > 140 


(5.14) 


Hence in the half plane Of 7 > —1/140, u 3 , v 3 and p 3 can become infinite 
at a finite time. This singularity is associated with either the unstable zero 
solution or an unstable limit cycle solution. 
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Some clarification of the regions (i) and (ii) defined by 

(i) 9o > goi , 0/l < -1/140 

(ii) the area bounded by 0/~f — —1/140 , the curve C and go = 0 

is required since here an unstable limit cycle solution exists but a finite time 
breakdown of the type discussed above cannot occur. The possibility of the 
existence of an alternative form for the algebraic singularity was investigated 
but none was found. However it can be shown that an alternative breakdown 
on an infinite time scale can occur. Here u 3 , v 3 and p are written 

{u 3 ,v 3 ) = {u 3 o(y),v 30 (y))e st + (u 3 i(y), v 3 i(y)) + 0{e~ st ) 
p- Poe st + Pi+0{e~ st ) 


(5.15a, 6) 


where s > 0 so that u 3 , v 3 and p go to infinity exponentially at the same rate. 
This exponential rate of growth satisfies the nonlinear eigenvalue problem: 



W 30 = v 30 = v' 30 = 0, y = 0, 1. 

7 [ uu 30 dy + 0{u3o(O) - £ 30 ( 1 )} = 0 
Jo 


(5.16) 


The latter condition ensures that the terms proportional to e 2st on the right 
hand side of (4.2) vanish identically. Otherwise p would go to infinity like 
e 2st and the structure (5.15) would not be possible. The eigenvalue s must of 
course be positive for (5.15) to be valid. It turns out that the eigenvalue s is 
closely related to the limit cycle solution and the existence of the finite time 
‘blow-up’ solution (5.7). 

Firstly it should be noted that when 9 / 7 tends to any value on the curve C 
of Figure 5 then {u 3 0 , v 3 o) tends to (po/p)(U, V) where {U,V) is the limit cycle 
solution determined by (5.2)-(5.4). The eigenvalue s of (5.16) goes through 
zero on C so solutions of the form (5.15) exist in a neighbourhood either 
above or below C. Secondly it can be shown that eigenvalues of (5.16) with 
s » 1 exist either immediately above or below the line 9/ 7 = —1/140. This 
result is not immediately obvious but follows by seeking a solution of (5.16) 
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with s » 1. The functions u 3 o, £30 for 0(1) values of y satisfy the same 
equations as the blow-up solution (u 30 , v 3 o) . Again it is found that sub-layers 
at y = 0, 1 are required in order to satisfy the no-slip condition there. The 
solution in these layers is in fact different from that for the ‘blow-up’ solution 
but the outcome is essentially the same. Thus it is found that eigenvalues of 
(5.16) exist for 9/ 7 sufficiently close to the line 0/7 = —1/140. This region 
is above or below 0/ 7 = —1/140 depending on the sign of the difference in 
the shear « 3 o(0) — u 30 (l) the lower and upper boundaries. The eigenvalues 

of (5.16) appear or disappear on the curve C and the line 0/ 7 = —1/140. A 
numerical investigation of (5.14) shows that for go < goi these eigenvalues are 
connected so that the exponential blow-up can occur in region (ii). Beyond 
Qo — 0oi the situation is more complex and the eigenvalues emanating from C 
and 9 / 7 = —1/140 are not always connected. In Figure 5 the regions where 
eigenvalues of (5.14) were found numerically to exist are shaded. 

Thus if it is not possible for u 3 , v 3 and p to become infinite in a finite 
time there exists the possibility of a breakdown in an infinite time. The only 
region of Figure 5 where neither kind of breakdown can occur is the area 
below C between go = 0 and go = <701- In this region the stable zero or limit 
cycle solutions are presumably set up as a result of any initial disturbances 

It should also be noted that another type of exponential blow-up is always 
available for go > goi- This type of breakdown has p = 0, (u 3 ,v 3 ) = 

eS *( w 30> v 3o) where s is the eigenvalue for a Gortler vortex. This type of 
breakdown is not always stable in the sense that a small perturbation of p 
will be forced by the exponentially growing Gortler vortex and will possibly 
grow. Thus it is likely that this kind of breakdown cannot occur in regions 
where either the finite time breakdown or the exponential breakdown (5.15) 
exists. 

There are regions in Figure 4 where the finite time blow up or either or 
both of the exponentially growing solutions can occur. It is not clear which 
of these would be the ultimate state for an arbitrary initial perturbation. 
However one would expect that a sufficiently complicated initial perturbation 
would, in general, lead to the finite time singularity. This was certainly found 
to be the case for the numerical integrations of (4.1), (4.5) to be reported 
later in this section. 

To determine the relevance of the different regions of Figure 5 to the in- 
teraction problem of §4 it is necessary to compute the constant 9 / 7 appearing 
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in (4.5). This is shown in Figure 4 as a function of fio for oblique waves in- 
clined at an 0(1) angle to the flow direction. It can be seen that 9/ 7 increases 
monotonically with j3o and 0/7 = —1/140 when f3o/ao ~ .8887. The finite 
time breakdown discussed above can therefore only occur for waves inclined 
at angles of more than 41.6 degrees to the flow direction. TS waves inclined 
at a lesser angle break down on an infinite time scale. 

Finally in this section the results of a numerical integration of (4.1), 
(4.5) are reported. These calculations were carried out in order to check on 
the breakdown structures described above. The integrations were performed 
by use of finite differences with the spatial derivatives on the left hand side 
of (4.1) evaluated implicitly. The terms on the right hand sides of (4.1) and 
(4.5) were treated explicitly. This scheme is stable for time step At, the same 
order of magnitude as the step length Ay. 

Firstly consider the situation when go < <701 so that the limit cycle is 
stable below C in Figure 4. Here we found that any nonzero initial disturbance 
decayed to zero if a < 0 and approached the limit cycle solution otherwise. 
For values of 9 / 7 > —1/140 we found that, dependent of a, a suffiently large 
initial disturbance terminated in a finite time singularity. If a was negative 
however a sufficiently small initial disturbance was found to decay to zero. 
For values of 9 / 7 between C and 9/ 7 = —1/140 a similar result was found 
except that the amplitudes were found to grow exponentially when breakdown 
occured. 

The above results are entirely consistent with our previous discussion. 
Similarly the numerical calculations for go > <701 were consistent with the 
prediction given above. Clearly a numerical investigation of (4.1) -(4.5) can 
describe only a finite number of initial disturbances but, in all the cases inves- 
tigated, we found that if both the algebraic and exponential ‘blow-ups’ were 
possible then except in special circumstances the algebraic type took place. 
The exception corresponded in all cases to an initial disturbance with the TS 
amplitude set equal to zero. 

6. Conclusion 

It has been shown in the previous sections that there are at least two 
types of nonlinear interaction possible where small amplitude TG longitudinal 
vortices and TS disturbances interact in curved channel flows. The weakest 
‘Type P interaction is governed by the coupled quadratic amplitude equation 
typical of resonant triad interactions. This interaction can only take place in 
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the presence of curvature and breaks down at a finite time with the amplitudes 
of both modes becoming infinite. 

At larger amplitudes the ‘Type IP interaction occurs and leads to the 
evolution equations (4.1), (4.5). The outcome of this interaction is again a 
singularity at a finite time if the wavefronts of the TS waves make an angle 
with the vortices which is greater than 41.6°. Perhaps the most significant 
property of this type of interaction is that it can also take place in a straight 
channel. In this situation the Gortler vortex would decay to zero in the 
absence of the oblique TS waves, thus the longitudinal vortex structure is then 
a direct consequence of these waves. Thus, in a straight channel, sufficiently 
large TS waves are capable of driving longitudinal vortices with precisely the 
same structure as those which would occur if the channel were curved. This 
suggests that if transition in channel flows is caused by interacting oblique 
waves, then waves making an angle of more than 41.6° with the flow direction 
will lead to a catastrophic change in the flow structure. 

In external flows it is well known that, in the initial stages of transition, 
lon gitudin al vortex structures can develop in the flowfield. (See, for example, 
Klebanoff, Tidstrom and Sargent (1962)). In addition, turbulent bursts in 
boundary layers are thought to originate from the spanwise locations where 
these vortex structures are most unstable. It follows from the present work 
that if boundary layer growth does not destroy the interaction mechanism 
discussed in this paper then there exists a mechanism for producing the lon- 
gitudinal vortices observed in the initial stages of of boundary layer transition. 
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Figure 1. 8 as a function of a calculated from the eigenrelation 
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Figure 2b. Sin <p as a function of a n for the first Gortler mode. 
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